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A far- infrared system measures the full complex Faraday angle, rotation as well as ellipticity, with an unprece- 
dented accuracy of 10/irad/T. The system operates on several far-infrared laser lines in the spectral range 
from 0.3 to 6 THz and produces results as a continuous function of temperature from 10 to 310K and applied 
fields between ± 8T. Materials successfully measured include GaAs 2-DEG heterostructures, various high 
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I. INTRODUCTION 

Hall effect measurements are a powerful tool routinely 
used in condensed matter physics. Extension of these 
measurements to optical frequencies provides dynamical 
characterization of the off-axis conductivity axy where 
optical propagation is in the z direction. This off-axis 
conductivity when combined with the longitudinal con- 
ductivity (Jxx from standard optical measurements pro- 
duces the full complex conductivity tensor along with the 
Hall angle defined as tan(6>i:/) = dxy/crxx- 

This additional information available from infrared 
(IR) Hall measurements has been particularly useful 
in investigating strongly interacting electron systems, 
materials often defying paradigms in condensed matter 
physics established in the twentieth century. Since inter- 
action effects can enter into axy differently than axxr^ 
the IR Hall effect provides new insights to strong correla- 
tion physics. Similar to axx^ Hall spectroscopy probes the 
energy scales of the system including energy gaps and car- 
rier scattering rates. It also obeys sum rules that permit a 
more global view of the effects of the interactions on the 
electronic structure of materials. ^^^^^ IR Hall measure- 
ments have led to several recent important findings: re- 
futing the spin-charge separation scenario for the anoma- 
lous behavior of the dc Hall effect in cuprates^*^*^ firmly 
establishing that the anomalous dc Hall effect in the para- 
magnetic state of n-type and p-type cuprates is caused by 
current vertex corrections produced by electron-electron 
scattering mediated by antiferromagnetic fiuctuations,^ — 
observing the first evidence of small Fermi pockets in un- 
derdoped cuprates, ^^ demonstrating that inelastic scat- 
tering enters into the Hall conductivity differently than 
the longitudinal conductivity,^^^ and characterizing the 
Fermi pockets formed by a spin density wave gap in un- 
derdoped Pr2-a;Cea;Cu04.^^^ 

Hall effect measurements at very low optical excitation 
energies in the THz spectral region probe intrinsic prop- 
erties associated with the Fermi surface (FS) which are 



sensitive to the FS topology as well as Fermi velocity and 
scattering rate anisotropics.^^ One of the advantages 
of the far-IR (FIR) Hall measurements over the mid-IR 
(MIR) Hall measurements^^ is that for many material 
systems, there are known intermediate energy scales be- 
tween 10 and 100 meV which make comparison between 
MIR Hall measurements difficult with dc transport and 
other low frequency probes. 

Important future applications include extending quan- 
tum Hall effect (QHE) measurements to the THz regime. 
The recent appearance of novel 2-D systems such as 
graphene and the conducting surface states on topolog- 
ical insulators (TI) has stimulated interest in the trans- 
port and magneto-transport in systems with "massless" 
Dirac electronic dispersion. Topological insulators are 
predicted to form a new 2-D quantum state of matter 
fundamentally related to the edge states in the QHE.^^- 
Although many exotic properties are predicted, of par- 
ticular interest are signatures at finite frequencies of a 
quantum Hall step axy = (l/2)e^ /h - even though there 
exist no edge states - never before observed. For that 
matter, the frequency dependence of the QHE has never 
been measured in any material at optical frequencies, a 
topic we are currently pursuing. In particular, the Fara- 
day angle in both graphene and TIs is predicted to be 
quantized in units of the fine structure constant^ — 

The instrument to be described in the following pages 
is capable of measuring the full complex Faraday angle 
defined as: 



tan^F 



(1) 



where tyx and txx are the transmission amplitudes 
through the sample, which is situated in a magnetic field 
perpendicular to its surface and parallel to the propa- 
gation direction of the optical beam. The system oper- 
ates on various laser lines in the terahertz spectral region 
(0.3 to 6 THz) and produces data to an unprecedented 
accuracy of 10/irad/T and as a continuous function of 
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temperature from 10 to 310 K. Using the full complex 
Faraday angle, together with independent measurements 
of the longitudinal conductivity (Jxxi one may then cal- 
culate both the the complex Hall angle and (Jxy 

In what follows we shall first present an overview of 
the system followed by a detailed examination of several 
of the more critical components and design issues. Fi- 
nally, we discuss the system calibration and operation. 
Useful mathematical notation and models appear in the 
appendices. 



II. INSTRUMENT OVERVIEW 

Fig.[T] presents a schematic representation of the sys- 
tem. At the beginning of the optical path, a chopped, 
continuous wave CO2 laser pumps a single longitudinal 
mode, far-infrared (FIR), molecular vapor laser which is 
the source of the THz radiation. An aluminum parabolic 
mirror focuses the output of the FIR laser onto a pinhole 
aperture to remove spatial noise. A second aluminum 
parabolic mirror collimates the emerging beam and di- 
rects it into a beam steering assembly consisting of two 
fiat mirrors. Next, an aperture stop reduces the beam di- 
ameter to a size so that the radiation does not encounter 
the sides of the bore tube of the mechanical rotator. A 
thin mylar sheet beam splitter splits off a small portion of 
this reduced diameter beam for use as a power reference. 
The main portion of the reduced diameter beam passes 
through the mylar beam splitter and enters the bore of 
the high speed mechanical rotator. A waveplate mounted 
within the bore typically spinning at ~3000 RPM modu- 
lates the beam so as to introduce various harmonics of the 
rotation frequency into the polarization states. A third 
aluminum parabolic mirror focuses the beam onto the 
sample located in a ±8 Tesla superconducting Helmholtz 
magnet. The magnetic field is aligned parallel to the 
beam path. 

In the Faraday geometry, the portion of the beam 
transmitted by the sample passes through an analyzer 
polarizer and onto a fourth aluminum parabolic mirror, 
which re-collimates it. A fifth aluminum parabolic mir- 
ror focuses the beam into the output detector. Note that 
the beams leading to and from the magnet and the fourth 
and fifth mirrors are collimated. This allows increasing 
the separation between the magnet and the detector and 
laser source so as to reduce the effects of the strong mag- 
netic field upon these components. Also note that, except 
for the small amount of spatial noise introduced by the 
aperture stop prior to the mechanical rotator, the system 
operates essentially within a single spatial mode. This 
limits noise from modal interference whose frequencies 
resulting from vibration, thermal drifting, etc. would 
otherwise fall within the mechanical rotator frequency 
and its harmonics. 

An ideal rotating waveplate transmits linearly polar- 
ized light when either the ordinary or extra-ordinary axis 
aligns with the polarizer. This occurs four times per rev- 



split coil magnet absorbers 
sample 



kapton 
polarizer 




FIR laser 

Au/Si/quartz 
output coupler 



reentrant Cu bore tube 
^^^^___ kapton 

"^ A 

>1.5m away 

waveplate 
■ I ^^^^^ — rotator 
polarizer 
beam splitter 
stop aperture 



beam steer 
aperture 



FIG. 1. (Color online) Schematic of the complex Faraday 
angle measurement system operating on various laser lines in 
the THz region. 



olution of the waveplate giving rise to a fourth harmonic 
signal. With a moments thought, one will recognize that 
a sample inserted in the magnet so as to cause a real 
Faraday rotation Re(^F) will change the relative phase 
of this fourth harmonic. In fact, the fourth harmonic 
phase change is directly proportional to the Faraday ro- 
tation. An angular position sensor on the mechanical 
rotator allows determination of this phase change. 

Considering no sample in place, the polarization state 
oscillates between pure left and right handed circular po- 
larized light twice upon one full rotation of an ideal wave- 
plate. Upon introduction of a sample, one can see quali- 
tatively that the second harmonic signal is then sensitive 
to the induced ellipticity or circular dichroism. The ap- 
pendix presents, in a very manageable notation, a com- 
plete derivation of all the various signals one might expect 
for typical samples. Corrections to accommodate errors 
and losses in the various system components is accom- 
plished via an in situ calibration technique discussed in 
Sec. UVCl 

In general, a sample induces polarization rotation, 
characterized by Re(^F), as well as ellipticity, charac- 
terized by Ivh^Of)' The resulting signals recorded at the 
output detector include dc level Rdc of the transmittance 
through the sample measured at the chop frequency uOchop 
of the CO2 pump laser and Ai^2oj, ^o,2a;, ^i,4u;, and Ao,4u;, 
which are the second and fourth harmonic signals on the 
output detector that are in-phase (i) and out-of-phase (o) 
with the original harmonic signal that results when no 
sample is present or when the magnetic field is zero. By 
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taking appropriate ratios of ^^,20;, ^o,2u;, ^i,4u;, ^o,4u;, 
and Rdc^ one may extract the complex Faraday angle 
via Eq. lA2[ Note that similar to heterodyne detection 
schemes, for small angles the desired quantities are pro- 
portional to the phase and therefore the transverse txy 
amplitude rather than being proportional to the power 
which is not so accurately measured for small values. 



III. CRITICAL COMPONENTS AND DESIGN ISSUES 



of Faraday experiments. 

In some very specific circumstances, etalons within 
the sample can be used to emphasize particular mate- 
rial properties. For example, appropriate engineering of 
the etalon in TIs allows isolation of the surface state re- 
sponse producing a Faraday signal which exactly equals 
the fine structure constant. ^^ However, in most circum- 
stances, elimination of the etalon is desirable to simplify 
the optical system and subsequent analysis. 



A. Etalons and Multiple Reflections 



B. Sample design and mounting considerations 



Multiple refiections within optical elements and be- 
tween optical elements can result in spurious signals and 
must be minimized. Consider, for example, that the 
Faraday effect accumulates regardless of the direction of 
the light beam. Whereas a single pass through a sample 
or optical element may produce a rotation of 6>, a pass 
forward, then backward, then again forward will produce 
a rotation of 3^. This is a concern not only for the sam- 
ple but for any optic within the magnetic field, for which 
reason we use thin polyimide film, 25 micron, as the 
magnet vacuum windows. Even outside of the magnetic 
field multiple refiections can cause small variations over 
time due to very small motions or vibrations of the hard- 
ware like fiexing of the magnet cryostat or fiexing of the 
the aluminum optical table from cryogen boil-off, motion 
of the fioor, thermal drift, etc. 

To control these effects, optical elements may be: cho- 
sen to be much thinner than the wavelength; conditioned 
with an antirefiection (AR) coating; and/or set at an an- 
gle to the optical path. For example, a dielectric AR 
layer was utilized on an LaSrGa04 (LSGO) substrate to 
impedance match between vacuum and the FIR substrate 
index of n ~ 4.4. The matching layer was a piece of z- 
cut quartz chemically etched to a thickness equal to an 
odd multiple of a quarter wavelength and held against 
the substrate by means of a spring. The refractive index 
of quartz in the FIR is 2.12 which is close to the geomet- 
ric mean, 2.10, of vacuum and LSGO. High resolution 
FTIR transmission measurements, corresponding to var- 
ious laser frequencies^^ were used to confirm the etched 
thickness. 

Another type of broadband AR coatings that may be 
applied to sample substrates, waveplates, and some de- 
tector cold filters is a thin metallic film.^^ The Fres- 
nel reflection coeflicient for an interface between vacuum 
and a thin metallic film on a dielectric of index n is 
r = {n — l — Zo/R[j)/{n-\-l-\-Zo/RQ) which is zero under 
the condition that the sheet resistance Rc\ = Zo/{n — 1) 
where Zq = 377 Q is the impedance of free space, i^n 
is approximately constant provided the scattering rate is 
substantially greater than the optical frequency, a condi- 
tion met in the FIR for nichrome (NiCr) films. Such AR 
coatings on quartz, silicon, gallium arsenide (GaAs), and 
LSGO were characterized via high resolution FTIR trans- 
mission measurements and used successfully in a number 



Ideally the sample should be made as large as possi- 
ble and the laser spot as small as possible. However, 
diffraction limits the minimum spot size at the sample 
to 1.22Xf/D. For our magnet cryostat, f/D is 3.5 so 
the spot size for 20 cm~^ radiation is 2 mm. Samples 
much smaller than the spot size can be located within an 
aperture to reduce the background. However, an aper- 
ture introduced at this point in the beam path will cause 
spurious polarization changes at its edges that can mask 
the Faraday signals. Constructing the aperture from an 
absorbing material like graphite reduces such effects. 

Many of the samples studied were highly refiective and 
poorly transmissive. By strategically placing graphite 
absorbers on the re-entrant copper bore tubes (Fig.[T]) 
and sample stick (Fig. [2]), and mounting the sample on 
an aperture, the effects of beam scatter were measured 
tobe~5xl0~^ which is negligible in nearly all circum- 
stances. Imaging the sample onto a spatial filter placed 
at the location of the detector has been used effectively to 
further attenuate scattering effects. Note that an aper- 
ture in this location will not introduce any spurious signal 
by modifying the polarization since the analyzer polarizer 
has already decoded the polarization information. 

For the simple case of thin film samples, the Faraday 
angle is related to the Hall angle as shown in Eq. lA15l bv a 
factor involving the sheet resistance l/{axxd) acquired by 
an independent FTIR spectroscopic measurement. This 
value may be engineered to be close to 1 by choosing 
the thickness of the film. However, as shown in Eq. lA15| 
the sheet resistance should be as small as possible while 
still allowing adequate transmission. As a gauge, sam- 
ples which attenuated the throughput intensity by four 
orders of magnitude were successfully measured (like the 
Bi2Sr2CaCu208+a^ sample deep in the superconducting 
state as shown in Fig. [5]).^ The signal at this level is de- 
tector noise limited. 

The sample mounting system must maintain the sam- 
ple position as well as control the sample temperature. 
Its central component is a sample stick, which was mod- 
ified from a design originally made by Oxford Instru- 
ments. The cold finger shown in Fig. [2] is cooled by a 
continuous fiow of liquid helium. The sample rests on a 
1/16" thick copper plate which has a graphite aperture 
smaller than the sample. A 100 1^, 1/4 W metal film 
resistor epoxied into a slot in the sample plate serves as 
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FIG. 2. (Color online) Photographs of the sample stick cold 
finger depict two views. Total length shown is ^2.5 inches. 
The top sample is peeled single crystal Bi2Sr2CaCu2 08+x (Bi- 
2212) on an AR coated quartz substrate. The bottom sample 
is a NiCr coated GaAs 2-DEG heterostructure with an AR 
coating. 



a heater. To restrict the motion of the sample during 
temperature and magnetic field sweeps, a 3" long thin- 
walled stainless-steel tube is silver soldered to the end of 
the cold finger. This tube fits into a teflon insert in an 
aluminum chuck which we added to the cryostat housing. 
The chuck is actuated from outside the housing. A cop- 
per wire serves as a thermal link whose size is adjusted to 
allow sufficient cooling of the sample while throttling heat 
flow to the sample stick from the sample heater minimiz- 
ing thermal expansion effects. The temperature can be 
continuously varied from 10 to 310 K. Thermal stresses 
on the sample are minimized while maintaining adequate 
thermal coupling to the copper plate by mounting the 
sample with a spring clip on one end and silver epoxy on 
the other. A Lakeshore 340 temperature controller mon- 
itors and regulates the temperature of the sample and 
simultaneously monitors the temperature of the sample- 
stick. 



C. Magnet 

Our magnet is an Oxford Instruments Helmholtz or 
split coil superconducting magnet that can scan mag- 
netic field between ±8 T in approximately 12 minutes. 
The magnet is maintained in a liquid helium bath. The 
vacuum windows are 0.001" thick polyimide which are 
slightly permeable to helium, allowing a small amount of 
exchange gas to enter the vacuum chamber. This causes 
cooling of the polyimide windows resulting in water con- 
densation. This is mitigated by blowing room tempera- 
ture air constantly across the windows. The entire mag- 
net dewar rests on four 1/4-20 screws to allow height 
and tilt adjustment for aligning to the optical beam. An 
Oxford model PS 120- 10 power supply controls the super- 
conducting magnet current. 



D. Laser system 

As mentioned earlier, the laser system is a molecular 
vapor FIR laser pumped by a continuous wave CO2 laser. 
The CO2 laser is a model PL5 from Edinburgh Instru- 
ments Ltd. and is tunable from 9 to 11 /im. Its output is 
a 7.5mm diameter beam with maximum power of SOW. 
The output is focused by an AR-coated ZnSe lens into 
the cavity of the FIR laser through a small hole in one 
of the FIR cavity mirrors. 

The FIR laser is a home-built dielectric waveguide 
laser— ii^ with two adjustable flat mirrors on either end. 
The input coupler is a 2 inch diameter molybdenum flat 
mirror with a 1.5 mm aperture purchased from SPAWR 
Industries. The input coupler accommodates a ZnSe 
Brewster window which minimizes destabilizing back- 
reflection feedback into the CO2 laser. The output cou- 
pler is a quartz substrate with a silicon dielectric coating 
grown to a thickness matching a 1/4 wave at 10 jam. to 
provide a high reflectivity for the pump beam. A gold 
annulus sputtered on the silicon surface defines a 1cm 
diameter FIR output aperture. The total area of the out- 
put coupler allows the CO2 beam to be maximally main- 
tained within the FIR laser cavity while concurrently fil- 
tering out MIR frequencies from optical elements down- 
stream. The large aperture minimizes diffraction effects. 
A 1.3 m fused silica tube forms the FIR waveguide and 
maintains its length. Its low thermal expansion coeffi- 
cient 7.5 X 10~^ /K helps to maintain stability. Diffrac- 
tion losses strongly attenuate higher order modes ^^ while 
lower order modes are only slightly attenuated due to 
the grazing incidence with the waveguide. ^^^—i^Sr— The 
power output depends strongly on the chosen FIR fre- 
quency, but is typically between 1 and 50mW. The axial 
mode spacing is Auax = c/(2L) ~ 150 MHz, the Doppler 
broadening is ~ 5 MHz for methanol, and pressure broad- 
ening is typically ~ 40MHz/torr ^ 4 MHz under the 
operating pressure of lOOmTorr. Typical gasses used 
are CD3OD (24.6cm-i), CH3OH (84.7cm-i), CH2F2 
(42.3 cm-^), and CH3OD (175.4 cm"^). A continuous 
flow at typically ~100 mTorr pressures is maintained, and 
a gas handling system with a cold trap is used to recap- 
ture the gasses. 

The FIR laser lines are scanned for the presence of 
multiple resonances using a step scan fourier transform 
polarizing spectrometer. The cyclotron frequency of a 
calibrated GaAs 2-DEG provides a convenient in situ tool 
to ascertain the FIR laser wavelength (see Sec. lIVD]) . 



E. Mechanical Rotator 

The high speed mechanical rotator spins an optic of up 
to 1" diameter at an adjustable speed up to 3000 rpm re- 
sulting in a rotator fundamental frequency of 50 Hz. The 
optic is located within a nominal 1" diameter threaded 
bore tube. Two threaded rings sandwich the optic be- 
tween two o-rings that provide cushioning to distribute 
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stress which would otherwise introduce non-uniformities 
into the optic. The bore tube is pressed into two ABEC- 
7 angular contact ball bearings which are themselves 
pressed into the aluminum rotator housing. An Animat- 
ics model SM2315D servomotor with programmable PID 
feedback drives the bore tube by means of a kevlar rein- 
forced polyurethane timing belt which engages a pulley 
pressed onto the bore tube end. An intermediate belt 
tensioner reduces belt slap and jitter. Black paint on a 
portion of the bore tube pulley provides a target for an 
optical sensor which supplies the electronic phase refer- 
ence for the lockin amplifiers. The phase noise as mea- 
sured on the lockins is within 1 : 10^ (at 50 Hz using a 
200msec time constant). 

An array of various thicknesses of x-cut quartz first- 
order waveplates are chosen to cover the entire spectral 
region from 3 to 240 cm~^. Whenever the thickness of 
the waveplate is not optimal resulting in a retardance of 
A/3 ^ 7r/2, signal is sacrificed (see Eq. lA2p . The lim- 
iting upper frequency when using quartz waveplates is 
- 240 cm-^ 



F. Detectors, mirrors, polarizers, and beam splitter 
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FIG. 3. (Color online) The real and imaginary part of the 
Faraday angle is measured as a function of angle for a quartz 
waveplate and fit using the waveplate parameters S/S and ^F. 
The imaginary part is calculated using the fit parameters and 
overlayed on top of the data. 



The detectors are silicon composite bolometers pur- 
chased from Infrared Laboratories. The primary detector 
is cooled to 2K and the reference detector to 4K. The 
77K and helium cold stage filters are quartz, the vacuum 
window is wedged white polyimide, and a Winston cone 
collects radiation up to an f/# of 3.5. The characteristic 
frequency response of the primary detector is ^500 Hz 
with a noise-equivalent-power of 2 x 10~^^ W/Hz^/^ at 
100 Hz and a responsivity of 2 x 10^ V/W. 

All 90 degree off-axis parabolic mirrors are 4" alu- 
minum substrates with MgF2 coatings purchased from 
Janos technology. Flat mirrors are bare aluminum evap- 
orated on optically flat glass. 

Polarizers are evaporated Au on 6 /am. thick mylar 
with 4 jam pitch and 2 jam. linewidth purchased from Sci- 
enceTech. These are mounted taut between aluminum 
rings. Their measured extinction coefficient is 10 ~^ at 
84.7 cm-^ 

The reference beam splitter is a 0.00048" thick mylar 
sheet that is mounted taut between 2 flanges with o-rings. 
The reflectance is ~10% at 84.7 cm~^ and is slightly fre- 
quency dependent. 



G. Electronics 

The bolometer detectors are equipped with pre- 
amplifiers that allow a high impedance (^ 50 1^) output 
at gain settings of 1, 200, or 1000. 

The output signal of the transmission bolometer de- 
tector is fed into an EG&G 5113 low- noise pre-amplifier 
which is then fed into a 7280 EG&G lock-in amplifier. 
A homemade optical encoder reads the frequency of the 



rotator bore and generates TTL pulses to the lock-in am- 
plifier. The 7280 decomposes the signal into in- and out- 
of phase components of the second and fourth harmon- 
ics. An EG&G 7260 and a EG&G 7265 lock-in amplifier 
measure the magnitude of the chopped signal from the 
reference and transmission detectors. 

The detectors are nonlinear causing sums and differ- 
ences of frequencies of various noises and signals to ap- 
pear in the spectrum. A spectrum analyzer enables 
choosing cOchop and u (typically 311 and 49 Hz, respec- 
tively) such that the frequencies of interest do not inad- 
vertently overlap with other harmonics or their sums and 
differences. 

An in-house Labview program simultaneously controls 
and acquires data from the magnet power supply, tem- 
perature controller, EG&G 5113 preamplifier, and three 
lock-in amplifiers. 



IV. INSTRUMENT CALIBRATION AND OPERATION 

A. Electronics 

1. Lock-in amplifiers: harmonic analysis and setting the 
phase 

The complete detector signal is given by Eq. lAll and 
[A2l The EG&G lock-in amplifiers measure the RMS- 
amplitude of all signals. Optically chopping produces a 
square wave whose first harmonic amplitude is 0.637 com- 
pared to the square wave amplitude. The measured Rdc 
value needs to be multiplied by a factor of 2 since the ac- 
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tual dc level is the peak-to-peak value of the square wave- 
form. In this case, the magnitude of the chopped signal 
measured by the lock-in amplifier needs to be multiplied 
by 2 /mi to be used in Eq.[A2]as Rdc 

The phase of the Auo lock-in is set to zero in zero mag- 
netic field. However, inspection of Eq. [A2l shows that the 
2uj phase can not accurately be set under the same con- 
ditions since both the in- and out- of phase signals are 
close to zero (since AF ^ 0). Conceptually, the 2uo phase 
can only be set when there exists a substantive Im(^i?) 
signal such that Ao^2uj » ^i,2a;- Since AF is small, max- 
imizing Im(6>i?) by adjusting the phase of the 2uj signal in 
magnetic field is, in principle, a good method for setting 
the phase. In practice, a single lock-in is used to measure 
both the 2uj and Auj signals. The phase of the second har- 
monic is set equal to the phase of the fourth harmonic. In 
software, the phase of the 2uj channel is then rotated by 
A(j)2u; to maximize the Im(6>F) signal. A^2u; is found to be 
the same for a variety of samples and corresponds to the 
phase shifts observed when characterizing the frequency 
dependence (where phase information was acquired) of 
the detector between the 2uj and Auj frequencies. 



2. Detector frequency roll-off 

The relative phase and amplitude response of the de- 
tectors and preamplifiers as a function of modulation fre- 
quency are experimentally characterized. One method 
involves double chopping an FIR beam whereby one 
chopper continuously sweeps frequency while the other 
remains at fixed frequency. The amplitude response is 
used to correct the detector signal between the three rel- 
evant modulation frequencies: 2a;, Auj^ and uOchop- 



B. A/3 and AF from FTIR spectroscopy 

FTIR spectroscopy was performed on waveplates in- 
serted at various angles 9 between two parallel polarizers. 
The transmitted intensity (relative to the input intensity) 
is given by the following: 

+ 2e-(^^+^2) cos(A/3) sin2((9) cos'^{0) 

where 71 and 72 are the absorption along the e- and o- 
axis of the waveplate, and A/3 is the retardance. Taking 
ratios of I(^) for various angles 6>, AF = 72 — 71 and A/3 
as a function of frequency can be derived. 



C. A/3 and AF in situ calibration: polarizer as sample 

The preferred method to calibrate waveplates is in situ 
in which the unknown waveplate is inserted into the ro- 
tator. Adjusting a polarizer to various angles a inserted 
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FIG. 4. (Color online) The Faraday angle measured at 3 
frequencies as a function of magnetic field for an AR-coated 
GaAs 2-DEG heterostructure with a 25 Q/B layer of NiCr 
at 80 K. The real part of the Faraday angle is fit using 
the mobility (/i), electron density (n), and NiCr electron 
scattering rate {'jNiCr) as parameters. For all frequencies, 
n = 1.78 X lO^^m"^ and jNiCr = 5100 cm~^. The mobility 
was found to be 14, 12, and 8.4 m'^/Vsec for 24.6, 42.6, and 
84.7 cm~^ , respectively. The GaAs electron mass was taken 
to be 0.07 X me. 



immediately downstream from the rotator yields a de- 
tector signal given by Eq. lA4[ After taking ratios of the 
various amplitudes, there are a total of four equations. If 
we treat as unknowns the four quantities A/3, AF, and 
the frequency roll-off of the detector between 2uj and Auj 
and between 2uj and uOchopi then simultaneously solving 
for all four quantities gives an in situ calibration of the 
waveplate parameters and frequency roll-off of the de- 
tector. Small polarization errors introduced from imper- 
fections associated with the waveplate (for example, an 
imperfect AR coating) and/or induced effects from the 
kapton vacuum windows are corrected with this method 
by giving an effective A/3 and AF which is always within 
~ 5% of independent characterizations of the waveplate. 
The detector frequency roll-off never requires correction 
of more than small fractions of a percent. An analysis 
justifying this method of error correcting is detailed in 
Ref.|2a. 



The real part of the Faraday angle measured by rotat- 
ing the polarizer to small angles (±5 degrees) by Eq. lA3l 
exactly equals a (given by Eq. lA4p thus providing a very 
good absolute calibration of Re(6>i7). 
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D. Calibration of lm(^i?): inserting a waveplate or 2-DEG 
as a reference sample 



By inserting a waveplate immediately downstream 
from the rotator at various angles a while measuring 
Re(^F) and Im(^i?) given by Eq. lA7| it is possible to cal- 
ibrate the imaginary part of the Faraday angle with re- 
spect to the real part. The inserted waveplate parameters 
A/3' and AF' are obtained by fitting the real part of the 
Faraday angle. The fit parameters are then used to calcu- 
late the imaginary part, and overlayed on the imaginary 
Faraday angle data. Results for the 84.7 cm~^ data using 
a 0.3 mm thick waveplate is shown in Fig.[3l Waveplates 
measured at other frequencies were checked in the same 
manner. The fit parameters of the waveplate agree with 
other independent waveplate characterizations (Sec. lIVBl 
and HVC]) . 

A GaAs 2-DEG sample is used for calibration of the 
Im(^i?) as well. Im(^i?) and HeiOp) have a well defined 
relationship as shown in Eq. lA12l allowing calibration of 
Im(6>F) with respect to Re(l9F). A GaAs 2-DEG het- 
erostructure with an AR coating and a 25 l^/D layer of 
NiCr on the opposite side of the wafer was deposited in 
order to attenuate the transmission level to ~ 1%. The 
transmission level and magnitude of the Faraday signals 
are similar to many of our experiments. The mobility and 
the density of electrons are used as fit parameters. The 
small Faraday response of the NiCr layer is accounted for 
with a simple Drude model where the electronic scatter- 
ing rate is a fit parameter. 

This sample was used to check the imaginary part of 
the Faraday angle with respect to the real part at three 
separate frequencies: 24.6, 42.3, and 84.7 cm~^ . The real 
part was fit, and the imaginary part was then calculated 
using the fit parameters from the real part. The results 
are shown in Fig.|4j 

The same electron density and NiCr scattering rate 
were used for all three frequencies, namely 1.78 x 
10^^771"^ and 5100c7Ti~^. The mobility is a function of 
frequency and found to be 14, 12, and 8.4 m^/ V sec 
at 24.6, 42.3, and 84.7 cm~^, respectively. The temper- 
ature was held at 80 K for all data. Since the mobility 
is inversely proportional to the effective mass as well as 
the scattering rate, the data suggest that either the mass 
or the hall scattering rate is enhanced as the frequency 
is increased. Since the resonant frequency is not signif- 
icantly shifted away from the calculated cyclotron fre- 
quency (oc 1/m), the frequency dependence is associated 
with the scattering rate. 

All of the data show some small asymmetry about the 
resonance in the imaginary part of the Faraday angle. 
This might be a result of some small amount of mix- 
ing between the real and imaginary parts of the Faraday 
angle. Imperfect sample AR coatings will cause etalons 
leading to some small amount of mixing. 



24.6 cm-^ 


42.3 cm-^ 


84.7 cm-^ 


Re(6'F) Im(6'F) 


Re(6'F) Im(6'F) 


Re(6'F) Im(6'F) 


2.031 
1.937 
2.057 
1.982 


0.374 
0.409 
0.394 
0.373 


1.992 
2.074 
1.992 
2.035 
2.051 


0.457 
0.450 
0.457 
0.446 
0.443 


1.865 
1.886 
1.957 
1.911 


1.081 
1.105 
1.067 
1.097 


2.002 


0.388 


2.029 


0.451 


1.905 


1.088 



TABLE I. Slopes of individual magnetic field sweeps in 
mrad/T performed on Bi-2212 at 100 K like those depicted 
in Fig.[5Ua,b). The last row is the column average. 



V. INSTRUMENT PERFORMANCE: MEASUREMENTS 
ON BI-2212 

To illustrate the performance of the instrument, we 
show reported results^ on optimally doped single crystal 
Bi2Sr2CaCu208+a, (Bi-2212) samples 100 nm thick with 
an area defined by an aperture 2.5 mm in diameter as 
shown in Fig.[2l The thin cuprate superconducting crys- 
tal was mounted to an AR coated quartz substrate. The 
longitudinal conductivity (Jxx was independently mea- 
sured via FTIR spectroscopy. 

Magnetic field sweeps to ±8T were performed at 
100 K. Fig.[5fa,b) shows the Faraday angle at 100 K as 
a function of magnetic field for 10.5 meV (84.7 cm~^) 
radiation. A compilation of the slopes associated with 
individual magnetic field scans of both the real and imag- 
inary part of the Faraday angle for all frequencies is pre- 
sented in Table [D 

Temperature sweeps from 100 to 300 K and from 300 to 
25 K are shown in Fig.[5fc,d). The transmission changes 
by two orders of magnitude between high and low tem- 
perature since the sample becomes superconducting. Ev- 
ident in the traces is the growing level of detector noise 
as temperature is decreased. There exists a tempera- 
ture dependent background produced by sample move- 
ment. Subtracting scans in ±8T and dividing by 2 nulls 
the background signal. The resulting Faraday angle is 
shown in Fig.[5fe,f) where multiple scans were averaged 
and then smoothed with a 3 K moving average. Magnetic 
field sweeps taken at various discrete temperatures give 
the same value of the Faraday angle as the temperature 
sweep data. 

The current limitation of the signal-to- noise level is im- 
posed by drift produced by sweeping temperature and to 
a lesser extent from sweeping magnetic field. The current 
noise from these effects is O.lmrad, or 10/irad/T. 

Variations of the system are currently underway, dc 
current measurements which depend upon the chirality 
of the incident light are planned which utilize the same 
optical polarization modulation technique. ^^ A Kerr ge- 
ometry is being configured since the technique is sensitive 
to a single surface, particularly important in measuring 
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FIG. 5. (Color online) (a,b) Faraday angle measured at 10.5 meV (84.7 cm~^) and 100 K versus B-field. Four scans are depicted. 
(c,d) Faraday angle measured at ±8T versus temperature. Ten scans are depicted. (e,f) Faraday angle versus temperature 
found by subtracting scans in ±8 T, averaging multiple scans together, and performing a 3 K rolling average with temperature. 



topological insulators. 
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Appendix A: Mathematical formalism 
1. General Notation 

The entire optical system can be described as a series 
of transfer functions which operate on the initial state of 
an incident electro-magnetic plane wave. 

Dirac notation is used to efficiently keep track of 
the polarization state as is described in Ref. 25 and 
Ref.yjJ. Under our convention whereby the radiation 
propagation is in the positive z direction, the elec- 
tric polarization vector \E) can be represented in ei- 
ther a linear (L) or circular (C) polarization basis as: 



{L\E) = 



and {C\E) 



where E± = E^x ^ iEyy. E_^ then represents 
the circular polarization with the electric field rotat- 
ing in a positive direction about the z axis in the 



lab frame. The main types of optical components 
of our system are: polarizer (P), waveplate (WP), 
and sample (S), which may be represented as follows: 

{L\P\L) = ^^ ^' 



{L\WP\L) = 



,*/3i 















xy 



{L\S\L) = 

^yx ^yy ^ 

Ap = Pi — 132 and AF = Fi — F2 are the differential 
phase shift and absorption between linearly polarized 
light along the ordinary and extra-ordinary optical axis of 
the waveplate. For an ideal quarter waveplate A/3 = ^ 
and AF = 0. Actual values are experimentally deter- 
mined by an in situ calibration technique (Sec. lIVC]) and 
verified by FTIR measurements (Sec. lIVB]) . txx^ txy^ tyx, 
and tyy are the transmission amplitudes for an arbitrary 
sample. Most of the measured samples are rotationally 
invariant meaning that txx = tyy and txy = 



-tyx- In this 



{L\So\L) 



^XX 

-t 



^xy 



xy 



f^x 



The rotation of a component by an angle (j) may be 
represented by an arbitrary angle (j) clockwise by the fol- 
lowing operator: 

^cos (j) — sin ( 
sin 6 cos ( 



(L|i?(0)|L) 



A rotating quarter-waveplate spinning at an angular 
frequency u about the positive z axis is represented by: 

{L\RWP{uj)\L) = {L\R{ujt)\L){L\WP\L){L\R{-ujt)\L) 
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2. Detector Signal Analysis 

a. Rotationally invariant Sample 

Using the above notation, the intensity at the output 
detector of the system shown in Fig.[Tl for a rotationahy 
invariant sample is: 

Id = \{L\P\L){L\So\L){L\RWPiuj)\L){L\P\L){L\Eo)f 

= Rdc + Ao^2uj sin 2ujt + Ai^2oj cos 2ujt (Al) 

+ Ao,4uj sin4a;t + Ai^^^ cosAut 

where, 

Rdc = e-^\U^\^-{{l-\td.neF\^) cos Af3^ 

(3+|tan(9Fp)coshAr) 

Ao,2uj = e-^|t^^p(/m(tanl9F)sinA/3+ 

i?e(tan6>F)sinhAr) 

A,,2c. = -e-^|t^^psinhAr (A2) 

Ao,4uj = e~^ \t XX f -Re{tSin Of) {cos AP - cosh AF) 
Aaoj = e~^\ta,a,f-{\tdinOFf - l)(cosA/3-coshAF) 

Ai,2oj, Ao,2oj, ^i,4u;, ^o,4u; are the second and fourth 
harmonic signals on the output detector that are in-phase 
(i) and out-of-phase (o) with the original harmonic signal 
that results when no sample is present or when the mag- 
netic field is zero. By taking various ratios of the mea- 
sured quantities in Eq. lA2| one may extract the quantities 
Re(tan6>i?) and Im(tan6>i?). 

In many circumstances, txy is zero in zero magnetic 
field and remains small even up to the maximum mag- 
netic field such that \txy\'^ <C \txx\'^ • In this small angle 
limit, the Faraday angle may be obtained from Eq. lA2l 



Reieph = 

Re{eF)2 = 
Im{eF)i = 

Im{eF)2 = 



lAo. 



Auj 



1 COS A/3 + 3 cosh AF Aq^aoj 

2 COS A/3 - cosh AF Rdc 
-i?e(i9F) sinhAF 

sin A/3 ^ 

1 cos A/3 + 3 cosh AF Ao^2uj 
4 sin A/3 R^ 

-Re{OF)smhAT 

sin A/3 
1 cosh AF — cos A/3 Ao^2uj 
4 sin A/3 Ai^^^j 



(A3) 

The preferred ratios for evaluating the Faraday angle 
involve the second and fourth harmonic amplitudes. The 
corrections involving the detector frequency dependence 
become more sever for the case of Rdc due to the com- 
paratively high chopping frequency. 



b. Polarizer or waveplate as a sample 

Inserting a polarizer at an angle a with respect to the 
analyzer in place of the sample, and replacing the rota- 
tionally invariant sample in Eq. lAll with {L\SP{a)\L) = 

{L\R{a)\L){L\P\L){L\R{-a)\L) yields the following de- 
tector signal: 

g-r 
Rdc = —r- cos^a (cos 2a cos A/3 + (2 + cos 2a) cosh AF) 

g-r 
Ao 2uj = -— cos^a sin2QfsinhAF 



M,2oj 



^o,^ 



■ cos^a (1 + cos 2a) sinh AF 

- cos^a (cos A/3 — cosh AF) sin 2a 



AiAoo = -r- cos^a (cos A/3 - cosh AF) cos 2a ( A4) 

Notice that Ao^4uj/Ai,4:Uj = tan2Qf. 

Inserting a waveplate in place of the sample defined by 
a retardance and differential absorption given by A/3' and 
AF', respectively, whose ordinary axis is at an angle a 
with respect to the analyzer yields the following transfer 
function 

{L\SWP{a)\L) = {L\R{a)\L){L\WP\L){L\R{-a)\L) 

(A5) 

Replacing {L\So\L) with {L\SWP{a)\L) in Eq,\M\ 
yields detector amplitudes in terms of A/3', AF', a, A/3, 
and AF. 

We make the following definition: 



{L\SWP{a)\L) 



f' f' 

^xx ^xy 

f f 

yx ""yy 



(A6) 



With this definition, t'^y = ty^, and if AF' is sufficiently 
small, then tyy ^ t^^. Under these circumstances, the 
Faraday angle may be derived from the sample waveplate 
parameters: 



^e(M=^(l-^^)- '''^ 



ItL.WA^ 



Im{eF) 



'XX \ ^H,4:Ujt 

Re{OF) sinh AT 



sin A/3 



r(i 



t 



xy' 



+' |2 



COS A/3 — cosh AF Aq, 



2uj 



sin A/3 



^i,4:UJ 



where Kyl^ /If^x? is given by 



(A7) 



(A8) 



I (cosh AF' - cos A/3') sin^ 2a 
e~^^' cos^ a + e^^' sin^ a -\- ^ cos A/3' sin^ 2a 



3. Faraday angle from a Drude thin film sample 



A Drude response in magnetic field is derived by using 
the classical force equation qE -\- ^v x B — ^v = m^ 
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together with the definition of current J = nqv = dE. 
The conductivity tensor is diagonal when represented in 
the circular polarization basis: 



Eq. lAlll and using Eq. lAlQI gives: 



tan^ 



= iZ' 



cr+ 



1 



{C\cj\C) 



m 'y — i{uj—0Jc) 









nq 



7 — i(u;+u;c) / 







t++t_ " 2 + Z'(cr+ + cr_) 

Converting the conductivity back into the linear polar- 
ization basis (Eq. lAll]) and simplifying gives: 



(A9) 

Since the matrix relating J and E is diagonal, the com- 
plex Fresnel transmission coefficient for a plane wave at 
normal incidence from vacuum through a thin conduct- 
ing film (A ^ film thickness) on a semi-infinite substrate 
of index n can be written as: 



tan6>F 



Z' 



a. 



xy 



'xy 



1 + Z'cFxx {Z'cJxxd) ^ ^l (Jx 



Substituting the definition of the Hall angle (Eq. lA14]) 
and rearranging yields the desired result: 

td^nOH = 1 + — 3 tan Op (A15) 

V Zo CFxx dJ 



t± = 



1 + n + Zocr±d 



(AlO) 



where Zq = 377 ft/n is the impedance of free space, and 
d is the thickness of the film. 

Utilizing the following conversions. 
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0'± = (Jxx ± ^ O-xy 
t'lt — ^xx ^ ^ ^xy 

and combining Eq JA9l and lAlQI gives: 



(All) 



tan^F 



hZUJc 



(7 — ioj){'y -\- K, — ioj) + uj'^ 



(A12) 



where k = ^Z', cv, = ^^, Z' = d Zo/{n + 1), Zq = 
377 r^/n, n is the index of the substrate, and d is the 
film thickness. /^ is a frequency which characterizes the 
radiation damping of the 2-DEG, a plasma mode which 

takpg into arrnnnt tVip rpHnrtinn nf tVip plprtrir fiplH insirlp 
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Two-dimensional electron gas heterostructures are 
usu ally characterized in terms of the electron mobility 
and electron density. The mobility can be defined as: 
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Inverting Eq JA13l and substituting 7 into Eq JA12l gives 
the Faraday angle in terms of mobility. 



4. Converting from Faraday angle to Hall angle in the 
thin film limit 



The Hall angle is defined as: 



tan Oh = -^^ 



(A14) 



Assuming a rotationally invariant sample in the thin- 
film limit, and converting Eq.[T]to the circular basis via 
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